Generalized switchable coupling for superconducting qubits using double resonance 
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We propose a method for switchable coupling between superconducting qubits using double res- 
onance. The inter-qubit coupling is achieved by applying near-resonant oscillating fields to the two 
qubits. The deviation from resonance relaxes the criterion of strong driving fields while still allow- 
ing for a fully entangling two-qubit gate. This method avoids some of the shortcomings of previous 
proposals for switchable coupling. We discuss the possible application of our proposal to a pair of 
inductively coupled flux qubits, and we consider the extension to phase qubits. 



I. INTRODUCTION 

Superconducting systems are among the most likely 
candidates for the implementation of quantum informa- 
tion processing applications In order to perform 
multi-qubit operations, one needs a reliable method for 
switchable coupling between the qubits, i.e. a coupling 
mechanism that can be easily turned on and off. Over the 
past few years, there have been several theoretical pro- 
posals to achieve that goal @, S, H S, [E 0, IE IE EE EL 
and initial experimental advances have been made [12 . 
EE H, EE Ea E3, EE E3, E3, |2j|- The early proposals 
involved performing fast changes in the qubit parameters 
and taking the qubits out of their so-called optimal points 
@, Q or using additional circuit elements 0, IE @ • Both 
approaches increase the complexity of the experimental 
setup and add noise to the system. Rigetti et al. pj pro- 
posed a switchable coupling mechanism that is turned on 
by applying resonant oscillating fields to the qubits and 
employing ideas inspired by the double-resonance physics 
known from nuclear magnetic resonance (NMR) [IE HH . 
In their proposal the qubits are kept at their optimal 
points throughout the experiment, neglecting the oscil- 
lating deviations caused by the driving fields. In spite 
of its appealing minimal reliance on additional circuit el- 
ements, that proposal requires the application of large 
driving fields. Other authors later proposed alternative 
mechanisms that avoided that limitation while still us- 
ing oscillating fields or oscillating circuit parameters to 
induce inter-qubit coupling [E IE EE EH • Those most 
recent proposals, however, suffer from some limitations 
of their own, e.g. not being usable at the optimal point 
[H or requiring additional circuit elements [E EE EH • 

Here we propose a generalized version of the double- 
resonance method where the constraint on the driving 
amplitudes is substantially milder than that required for 
the proposal of Ref. @- Our proposal provides an al- 
ternative to experimentalists when deciding what is the 



most suitable coupling mechanism to use in their exper- 
imental setup. 

It is worth noting from the outset that the term double 
resonance could be somewhat misleading in this context, 
since the mechanism discussed below requires only one 
resonance criterion, namely the one given in Eq. ([6]). 
However, we use it following similar mechanisms in the 
context of NMR 

The paper is organized as follows. In Sec. II we intro- 
duce the model system and review recent proposals for 
achieving switchable coupling. In Sec. Ill we derive our 
proposed coupling mechanism and consider some aspects 
of its operation. In Sec. IV we discuss the possible ap- 
plication of the proposal to realistic experimental setups 
that use inductively coupled flux qubits or capacitively 
coupled phase qubits. We conclude the discussion in Sec. 
V. 



II. MODEL SYSTEM AND PREVIOUS 
PROPOSALS 

We start by describing the system in general terms, and 
we defer the discussion of its physical implementation to 
Sec. IV. The system that we consider is composed of two 
qubits with fixed bias and interaction parameters. Oscil- 
lating external fields can then be applied to the system in 
order to perform the different gate operations. In other 
words, we consider the same system that was considered 
in Ref. 
given by 



7]. The effective Hamiltonian of the system is 



H = 
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where ujj is the energy splitting between the two states 
of the qubit labelled with the index j; Q 



and 



'rj 



are, respectively, the amplitude, frequency and phase of 
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the applied oscillating fields, A is the inter-qubit coupling 

(i) 

strength, and CTq are the Pauli matrices with a = x,y,z 
and j = 1,2. The eigenstates of a z are denoted by \g) 
and |e), with a z \g) = \g). Note that we shall set % = 1 
throughout this paper. 

In order for the qubits to be effectively decoupled in 
the absence of driving by the oscillating fields, we take 
A <C A, where A = ojj — o-> 2 , and we have assumed, 
with no loss of generality, that io\ > 0J2 and A > 0. 
Note that the absence of terms of the form &i is 
also crucial to ensure effective decoupling. Let us also 
take A <C w, where uj represents the typical size of the 
parameters ujj. Since we will generally assume driving 
amplitudes £!j comparable to A, the above condition will 
be crucial in neglecting the fast-rotating terms below, i.e. 
in making the rotating- wave approximation (RWA). 

Single-qubit operations can be performed straightfor- 
wardly by a combination of letting the qubit state evolve 
freely, i.e. with Qj — 0, and irradiating it at its reso- 
nance frequency, i.e. taking ufi? = ujj. Under the effect 
of resonant irradiation, Rabi oscillations in the state of 
the qubit occur with frequency fij. 

Although a clear review of previous proposals is not 
possible without a detailed discussion, we summarize the 
ideas of those proposals briefly here. The proposal of Ref. 
involves irradiating each of the two interacting qubits 
on resonance, i.e. taking ojj = a*,-, and relies on one 
manifestation of double resonance [22|, [23| • The idea of 
the double resonance in that case is that not only is each 
qubit driven resonantly, but also the sum of the Rabi fre- 
quencies of the two qubits matches the difference between 
their characteristic frequencies (i.e. O1 + O2 = A). After 
making two rotating-frame transformations and neglect- 
ing fast-rotating terms, i.e. performing two RWAs, one 
finds that the inter-qubit coupling term is no longer ef- 
fectively turned off (note that those transformations are 
essentially a special case of the ones we shall give in Sec. 
III). One thus achieves switchable coupling between the 
qubits. That proposal was criticized, however, for requir- 
ing such large Rabi frequencies. The proposal of Ref. @ 
uses an external field applied to one qubit at the sum 
of or difference between the characteristic frequencies of 
the two qubits in order to perform gate operations (e.g. 
u>^ = u)\ — u>2> ^2 = 0). However, since all the relevant 
matrix elements, e.g. (#<?|<3x |ee), with the eigenstates of 
the Hamiltonian in Eq. (TTJ) vanish, the proposed method 
would not drive the intended transitions. One therefore 
needs to use a somewhat modified Hamiltonian, e.g. one 
that contains an additional single-qubit static term with 
a a x operator . In practice, that means biasing one of the 
qubits away from its optimal point in the case of charge 
or flux qubits. Since optimal-point operation is highly 
desirable in order to minimize decoherence, an alterna- 
tive mechanism was proposed in Refs. [9l In those 
proposals an additional circuit element that can medi- 
ate coupling between the qubits is added to the circuit 
design. That addition effectively makes the parameter 
A in Eq. (fTJ) tunable, with its value depending on the 



bias parameters of the additional circuit element. One of 
those parameters is then modulated at a frequency that 
matches either the sum of or difference between the char- 
acteristic qubit frequencies. Clearly, since the driving 
term contains the operator a^&x, it can drive oscilla- 
tions in the transition \gg) <-> |ee) or \ge) <-> \eg} 7 even 
when both qubits are operated at their optimal points. 
As mentioned above, however, the use of additional cir- 
cuit elements is undesirable, because of the increased cir- 
cuit complexity and decoherence. 

In the next section, we shall derive our proposal to 
couple the qubits by applying two external fields close 
to resonance with the interacting pair of qubits such 
that neither qubit is driven resonantly, but the sum of 
the (nonresonant) Rabi frequencies satisfies the double- 
resonance condition. Therefore, in some sense we relax 
the requirement that the driving amplitudes must be as 
large as A/2 , as is the case in Ref. Q , and we make up for 
the resulting loss of frequency by adding the qubit-field 
frequency detuning to the double-resonance condition. 



III. THEORETICAL ANALYSIS 

We now turn to the main proposal of this paper, 
namely driving oscillations between the states \gg) and 
|ee) by employing double resonance with non-resonant 
oscillating fields. We take the Hamiltonian in Eq. ([1]) 
and transform it as follows: 



H' = S\{t)HSi(t) 



dS\ 



where 



S\ (t) = exp 



^ 2 



(2) 



(3) 



A solution of the Schrodinger equation id\^> (t)) / dt = 
H\V(t)) can then be expressed as Si(t)\ *'(*)), where |*') 
satisfies the equation id\^'(t))/dt = H'\$'(t)}. To sim- 
plify the following algebra, we take (fii = <f2 = 0. Ne- 
glecting terms that oscillate with frequency of the order 



of uij , we find that 



3=1 



2 x 



+ - ^0'i 1) 0'i 2) COS&J rf t + COS 5uJ l{ t 

+4 1) 4 2) sin 5uj rt t - &PaW sin<W rf i ) ,(4) 



uj^ . We now make 



where 8uoj — ujj —Wj, and <5w r f = 
a basis transformation in spin space from the operators a 
to the operators f such that the time-independent terms 



3 



in Eq. ([4]) are parallel to the new z-axis and the y-axis is 
not affected. Equation 0$ can then be re-expressed as: 

3=1 V 7 

+~ ^t^t^ 2 ) cos 6>i cos 02 cos 5o; r ft + r!jpt!jp cos <L; r ft 
+r^ 1 ^f^ cos 2 smdujrft - tjpr^ cos 9i sinferfi 
+A (5) 



where Qj = ySujj + £1 2 , the angles 6j are defined by the 

criterion tan0j = Qj/Sujj, and A contains terms in Eq. 
(H]) that were not written out explicitly in Eq. ([5]) because 
they will soon be neglected. We now take the frequencies 
to match the criterion wi + u)% = 8lo v {, or more explicitly 

+ ttf + \JSojI + 0| = A - Slui + 5uj 2i (6) 

where, as mentioned above, 5uj = uij — o^ f , and A = 
u>i — u>2- We also take the two terms on the left-hand 
side of Eq. © to be comparable to one another. Taking 
the above condition allows us to simplify H' with one 
more transformation. Using a similar procedure to that 
we used above for the first transformation, we now take 

S 2 =expLp^f(Jh\, (7) 

and after neglecting terms that oscillate with frequency 
of order A we find that 

H" = A(i - cos 00(1 + cos0 2 ) {fWf W - fWfi 2 )} . 

The reason why we can neglect the term A in the above 
transformation can be seen by observing that all the 
terms contained in A contain at least one f z operator, 
and they oscillate with frequency <W r f. Therefore, even 
after the transformation, those terms will still oscillate 
with frequencies that are of the order of A (and ampli- 
tudes smaller than A), meaning that their effects on the 
dynamics can be neglected in H" , whose typical energy 
scale is a fraction of A. 

Equations © and (jHJl form the basis for the coupling 
mechanism that we propose in this paper. The Hamilto- 
nian H" drives the transition \gg) <-> |ee) but does not 
affect the states \ge) and \eg) in the basis of the oper- 
ators f . Therefore, a single two-qubit gate that can be 
performed using the Hamiltonian H" and the set of all 
single-qubit transformations form a universal set of gates 
for quantum computing. Note that since the two-qubit 
gate is performed in the basis of the f matrices rather 
than the & matrices, one needs to include in the pulse 
sequence the appropriate single-qubit operations before 



and after the two-qubit gate. Note also that if we take 
the special case cos0i = cos 02 = 0, i.e. 8u>i = 5^2 = 0, 
we recover the corresponding case in the results of Rcf. 

0- 

A first look at Eq. ([5]) shows that one can achieve faster 
gate operation than in the special case cos 6\ — cos 02 = 
by choosing cos0i to be negative and cos0 2 to be pos- 
itive. In other words, instead of using the special case 
of resonant driving (5ui — 5u)2 — 0) one chooses 5ui to 
be negative and 5u>2 to be positive (i.e., lu\ ! > w\ and 
UJ2 < uj 2 ). However, inspection of Eq. ([6]) while not- 
ing that J Slu 2 + f2 2 < 1 8u>j | + £lj shows that one would 

then have to increase at least one of the frequencies ttj 
above the value A/2 in order to satisfy Eq. ^ with that 
choice of 8u>\ and 8u)2- Since we started with the moti- 
vation of finding an alternative double-resonance method 
that works with smaller values of Qj , we focus on the op- 
posite case, namely 8lu\ > and 8u> 2 < 0, and we accept 
the resulting reduction in gate operation speed. Starting 
from the special case 8u>\ = 8ll>2 — and moving in the 
direction given above, we find that both f2s can now be 
reduced below the value A/2 while satisfying Eq. ©. 

It is worth pausing here to comment on the higher- 
order effects that we have neglected in making the two 
RWAs. The second-order shifts that we have neglected in 
making our first RWA, i.e. the Bloch-Siegert shifts, are 
of order ri 2 /wj [24|. That energy scale is not obviously 
smaller than the inter-qubit coupling strength A. One 
might therefore suspect that those shifts will prohibit the 
performance of the proposed method. That is not the 
case, however, since those shifts only modify the values 
of the required driving frequencies and amplitudes, as we 
shall demonstrate with numerical simulations in Sec. IV. 
There we shall take the case where A 2 /wi = 2A, and we 
shall show that full oscillations between the states \gg) 
and |ee) can still be obtained when the shifts are properly 
taken into account. Other frequency shifts that result 
from our approximations, and possibly other experiment- 
specific shifts, also affect the required driving frequencies 
and amplitudes. We will not attempt to give analytic 
expressions for those shifts. However, we will take them 
into account by numerically scanning the driving ampli- 
tudes to achieve optimal gate operation. 

We now ask the question of how low can VLj be cho- 
sen to be. In principle, Eq. © can still be satisfied 
by taking fl,- to be very small and taking 8lo\ w A/4, 
8lo2 ~ — A/4. Note, however, that the frequency of gate 
operations is given by the coefficient in Eq. |jHJ), namely 
A(l — cos0i)(l + cos02)/16. That coefficient therefore de- 
termines the width of the resonance, or in other words, 
the error tolerance in driving amplitudes from the reso- 
nance criterion (Eq. [S]) [25( . One is therefore restricted to 
using values of 0i and 02 such that the above coefficient 
is larger than the accuracy of the available pulse gener- 
ators. Furthermore, taking the inverse of the frequency 
determines the period of oscillations in the doubly rotat- 
ing frame, or in other words, the time required to per- 
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form a two-qubit gate operation. Since the decoherence 
time sets an upper limit on how slowly one can perform 
the gate operations, that consideration provides another 
restriction on the allowed values of 0\ and 62- An exper- 
imentalist must therefore take the two above considera- 
tions into account, along with any restriction they have 
on the maximum usable driving amplitudes, in order to 
determine the window of parameters where the coupling 
mechanism can be realized. The parameters can then be 
fine-tuned within that window for optimal results. 

As an added perspective to help visualize the reso- 
nance criterion, we show in Fig. 1 the relevant energy- 
level structure. One can compare this figure to Fig. 2 in 
Ref. [3] • In that case, the on-resonance Rabi frequencies 
provide all of the energy splitting (i.e. Cj\ and C02) re- 
quired to satisfy the resonance criterion. In the present 
case, the energy levels involved in the frequency matching 
are already brought closer to each other by the facts that 
(1) the difference wj f — u>2 is smaller than the difference 
u)\ — L02 and (2) the detuning of each driving field from 
its corresponding qubit brings the relevant levels even 
closer to each other. It would appear from Fig. 1 that 
the resonance criterion can be satisfied with arbitrarily 
small driving amplitudes and the proper choice of ui\ { 
and u>\ ■ As was discussed above, however, the matrix el- 
ement (in the dressed-state picture) coupling the relevant 
energy levels becomes very small in that case, leading to 
the undesirable situation of high required accuracy in the 
driving fields and slow gate operation. 
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FIG. 1: (color online) The energy level diagrams of the two 
qubits in the dressed-state picture. The resonance criterion 
is satisfied when the smallest energy difference between two 
adjacent manifolds of qubit 1 states becomes equal to the 
largest energy difference between two adja cent manifolds of 
qubit 2 states. Note that Cbj = y^tj| 



We reiterate that care must be taken in using the term 
double resonance in describing the coupling mechanism 
discussed above. However, since it seems that the term is 
used to describe a number of distinct phenomena in NMR 



[23| , some of which bear resemblance to the one discussed 
here, we have followed that broad definition of the term. 
Note, in particular, that the mechanism discussed above 
requires only one resonance condition, namely the one 
given in Eq. ([6]) . Neither applied field has to be resonant 
with its corresponding qubit, provided that they are kept 
close enough to resonance that the two-qubit gate can be 
performed in reasonable time. 



IV. EXPERIMENTAL CONSIDERATIONS 




FIG. 2: Two inductively coupled flux qubits. The symbols 
x represent Josephson junctions. The static and oscillating 
externally applied magnetic fluxes, and & T f(t), are used 
to control the two qubits. The interaction is mediated by the 
mutual inductance M between the two qubit loops. 



In the above discussion, we have not specified what 
kind of qubits we consider. Our results therefore ap- 
ply to any kind of qubit where the effective Hamiltonian 
of Eq. ^ describes the two-qubit system. Because of 
its relevance to current experimental attempts to achieve 
switchable coupling between superconducting qubits, we 
now focus on the case of two inductively coupled flux 
qubits, as shown in Fig. 2 [2(3, H(|. Since the trunca- 
tion of the full Hamiltonian to the effective Hamiltonian 
of Eq. ([I]) has already been discussed by several authors 
(see e.g. Ref. [8j) and it is not central to our discussion, 
we do not include it here. 

In experiments on flux qubits, the individual qubits 
typically have uij (27r) x 5 GHz (note that the exact 
value is not completely controllable during fabrication, 
with the uncertainty reaching 0.5-1 GHz in some exper- 
iments) [2(| [H, [I?]] • The inter-qubit coupling strength 
A can be taken to be around (2ir) x 0.1 GHz. The 
highest achievable on-resonance Rabi frequencies f2j are 
in the range of several hundred MHz to 1 GHz (times 
2tt). The achievable Rabi frequencies are therefore large 
enough when compared with the naturally (i.e., uncon- 
trollably) occurring inter-qubit detuning A, suggesting 
that it might be possible to implement the proposal of 
Ref. with the above qubit design. However, additional 
difficulties that we have not discussed in Sec. Ill arise in 
different experimental setups. 

One experimental difficulty arises when A is 0.5-1 GHz 
[2^ |. In that case, the required Rabi frequencies are large 
enough to excite higher states outside the truncated qubit 
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basis, in addition to exciting other modes in the circuit. 
One would therefore ideally want to avoid using the high- 
est values oiilj cited above (~ 0.5 GHz). Taking interme- 
diate values of cos 9 between and 1, the required Rabi 
frequencies can be reduced substantially, and the two- 
qubit gate operation can still be performed in a time of 
the order of a few hundred nanoseconds. That time scale 
is smaller than the qubit decoherence times (typically 
1-3 /is), which means that a simple two-qubit quantum 
gate operation could be observable in the near future. 
Clearly, an increase in the decoherence times would be 
highly desirable in order to achieve longer sequences of 
gate operations. 
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FIG. 3: (color online) The occupation probabilities of the 
four eigenstates as functions of time. The blue (black), green 
(gray), cyan and yellow lines (the last two are essentially zero 
and barely visible) correspond, respectively, to the states |ee), 
\d e ) an d \ e g)- The initial state is \gg), u)i/2n=5 GHz, 
W2/27r=4 GHz, and \/2n=0.1 GHz. The driving frequencies 
and amplitudes include shifts caused by higher-order correc- 
tions. In both (a) and (b), uif includes the Bloch-Siegert 
shift fi|/4a> 3 -. In (a) 8\ — ir — 82 = 7r/3, and the fis (approx- 
imately 2-7TX 0.29 GHz) were shifted by 0.5% to correct for 
shifts in our second RWA. In (b) 81 = 11 — 82 = 7r/8, and the 
fis (approximately 2tvx 0.1 GHz) were shifted by 6.38% [28| . 

We have performed numerical simulations to show that 
the two-qubit gate can be performed for a wide range 
of values of 6% and O2 (note that smaller values of 9\ 
correspond to smaller driving amplitudes, and that we 
take 6 = 9\ = tt — 62). The simulations are performed 
by solving the time-dependent Schrodinger equation with 
the Hamiltonian of Eq. ([1]). We therefore make the two- 
level system approximation in describing each qubit. The 
results are shown in Fig. 3. If we take realistic experi- 
mental parameters and 9 — 7r/3, which corresponds to a 
reduction in the required driving amplitudes by a factor 
of about two, and we take the qubit to be initially in the 



state \gg), we can see that the occupation probability os- 
cillates between the states \gg) and |ee) with negligible 
errors and a very reasonable oscillation period (note that 
since we are considering a simple experiment designed 
to provide a proof-of-principle demonstration of switch- 
able coupling, errors of the order of 1% are negligible). 
In Fig. 3(b), we take the same experimental parameters, 
but we now take 9 = 7r/8, which corresponds to a re- 
duction in the required driving amplitudes by a factor of 
five. We can see that full oscillations can still be achieved 
when taking into account the shifts in the required driv- 
ing fields. However, the period of oscillations and the 
required accuracy in tuning the driving amplitude are 
now outside the experimentally desirable range. These 
results therefore agree with the statement made above 
that one should look for the ideal point of gate opera- 
tion, i.e. reduce the amplitudes of the driving fields just 
enough to reduce the errors caused by them to acceptable 
levels. 
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FIG. 4: (color online) Same as in Fig. 3, but including the 
effects of 100% crosstalk (the occupation probabilities of the 
states \ge) and \eg) are now more visible than in Fig. 3, 
but they are still small compared to those of the states \gg) 
and \ee)). In (a) 81 — 82 = tt/2, the co ri s were shifted by 
0.5%, and the fis (2ttx 0.49 GHz) were shifted by 2%. In (b) 
81 = ty — 82 = 7r/3, the o/ f s were shifted by 2.5%, and the 
fis (approximately 2ttx 0.29 GHz) do not include any shifts 
from the expressions of Sec. III. 

Another experimental issue that we have not addressed 
above arises in the case of crosstalk, i.e. when each 
qubit feels the microwave signal intended for the other 
qubit [20(. In other words, the Hamiltonian describ- 
ing the system includes additional terms of the form 

(3Q,j cos(o^ f £ + (pj)(Jx , where j ^ j', and the coefficient 
(3 quantifies the amount of crosstalk. If the amplitudes 
of the applied fields are small, a microwave signal that is 
resonant with one qubit will not affect the other qubit. 
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However, if the Rabi frequencies are comparable to the 
inter-qubit detuning, e.g. when f2 2 = (uii —oj^)/^ an( i 
/3 ~ 1, crosstalk cannot be neglected. In our method the 
ratio fi 2 /(wi — 0J2) is equal to sin 9/ (2 + cos 9), suggest- 
ing that the harmful effects of crosstalk could be reduced 
by decreasing 9. In fact, we have verified with numer- 
ical simulations that the errors caused by crosstalk are 
reduced by using our method, as shown in Fig. 4. Some 
of the shifts to the driving frequencies and amplitudes 
were determined manually by looking for optimal results. 
Note that the driving parameters corresponding to Fig. 
4(b) also drive oscillations between the states \eg) and 
\ge). However, combining the two driven transitions still 
describes effective coupling between the qubits. The pe- 
riod of oscillations in Fig. 4(b) is about 100 ns, suggest- 
ing that an experimental demonstration of the coupling 
should be possible even in the presence of 100% crosstalk. 

Finally, let us make a few remarks about the possible 
implementation of our method to capacitively coupled 
phase qubits [3, [l|[ . It is perhaps clearest to start by 
noting a point that is not directly related to the proce- 
dure of implementing our proposal: one of the main con- 
siderations in charge and flux qubits, namely the ques- 
tion of optimal-point operation, is rather irrelevant to 
the study of phase qubits, at least in the usual sense 
of using eigenstates with special symmetries to minimize 
decoherence. The phase qubit is simply a single Joseph- 
son junction controlled by a bias current. The static 
part of the bias current determines the qubit splittings 
ujj, whereas the amplitude of the oscillating part of the 
bias current determines the Rabi frequencies Qj [29|]. If 
one now takes two capacitively coupled phase qubits, 
one finds that the coupling term has the form ay 1 '' ] cffi 
[30l ] . If we now take the phases of the oscillating fields 
<Pi = f2 = t/2, we can follow the derivation of Sec. 
Ill and obtain the same results. In phase qubits the 
qubit splittings ujj are typically a few GHz (times 2tt), 
and unlike flux qubits those splittings can be tuned using 
the bias current during the experiment. Rabi frequencies 
can reach a few hundred MHz, and the coupling strength 
can be taken to be (2ir)x 0.1 GHz, giving essentially the 
same values for the parameters as discussed above for 



flux qubits. We finally note that the driving fields are 
supplied through the bias current rather than through 
external fields, which means that crosstalk is not a prob- 
lem with phase qubits. Realization of our proposal, or 
even that of Ref. [7|, should therefore be possible with 
capacitively coupled phase qubits. 



V. CONCLUSION 

We have derived a generalized double-resonance 
method for switchable coupling between qubits. The 
qubits are driven close to resonance such that the sum of 
their Rabi frequencies is equal to the difference between 
the frequencies of the driving fields. Our proposal with 
nonresonant driving of the qubits relaxes the constraint 
on the resonant-driving proposal, i.e. that of Ref. 0], 
requiring large driving amplitudes. We have compared 
the operation of resonant and nonresonant driving. Al- 
though our proposal can be applied to any kind of qubits, 
we have discussed in some detail its possible application 
to the special, but experimentally relevant, case of induc- 
tively coupled superconducting flux qubits. We have also 
considered the possible extension to the case of capaci- 
tively coupled phase qubits. 
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